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PROGRAM DESCRIPTION 
PURPOSE 

METHOD 



This program performs the analysis of a reinforced 
concrete column under axial loading and moment in 
two directions (two-way moment). The program 
will also handle the analysis of a column under axial 
Load and bending in one direction (one-way moment). 

The program goes through four general phases . The 
first is the location of the XX and YY axis (see fig.l). 
The second is the computation of moments of inertia 
about the XX and YY axis . The third is the calculation 
of stresses for the four corners of the column and for 
the steel bar under maximum tensile stress. Phase 
four is the calculation of the neutral axis location 
based upon the stresses computed in phase three. 
The initial assumption for the neutral axis may differ 
greatly from the location computed in phase four. 
The user therefore has the option of a new analysis 
using the neutral axis location just computed as the 
assumed neutral axis location. This may be continued 
as often as necessary until close agreement is secured 
between the assumed neutral axis location and the 
computed neutral axis location. 

Input data for the analysis of a column is contained 

in three cards (see Page one ) . 

Card one contains the following information: 

1. X and Y dimensions of the column (see fig.l)col. may 
be square or rectangular. 

2 . Data covering location of steel reinforcing bars . 
Cover distance and bar spacing in X and Y direct- 
ions are given here. Bar spacing is assumed to 
be the same along AB and DC edges (see fig.l). 
Bar spacing is assumed to be the same along AD 
and BC edges . One row of reinforcing is assumed 
along all four column faces (see fig.l). 

3 . The assumed location of the XX and YY axis . 
Card two contains the following data: 

1 . Dimensions defining the points where the assumed 
neutral axis crosses the column edges (see fig.l). 

2. The area of the reinforcing used. This is the 
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area of one bar. All reinforcing is assumed 
to be uniform in size. 

3 . The axial load and the bending moments 
about the locations assumed for the XX and 
YY axis . 

4 . The number of reinforcing bars in the X 
direction along column edges DC and AB 
and the number of reinforcing bars in Y 
direction along column edges AD and BC . 
Corner bars are included in both counts . In 
figure 1 the column shown has 12 bars in 
both the X and Y directions . 

Card three contains the following: 

1 . The ratio of flexural modulus of elasticity 
for reinforcing steel to flexural modulus of 
elasticity for concrete and the fore-mentioned 
ratio minus one. 

The effect of steel reinforcing bars upon the 
transformed section varies according to whether 
the bar in question is in tension or compression . 
In the first case the ratio used is N, in the second 
case N-l is used. To facilitate this determination , 
the dimensions RPX1 ,RPY1 ,RPX2 , and RPY2 are com- 
puted (see fig. 2). RPX2 and RPY2 specify the 
location of point P2 . P2 is the point where the 
assumed neutral axis crosses the center line of the 
steel reinforcing near point B . RPX1 and RPY1 specify 
the location of point PI . PI is the point where the 
assumed neutral axis crosses the center line of the 
steel reinforcing near point D. 

The program next performs the computation for find- 
ing the sum of the transformed area due to all rein- 
forcing bars and the area moments about AD and DC . 
This is an iterative type procedure, beginning with 
the bar at corner D. Each bar is handled separately. 

Row one reinforcing is considered first. The number 
of reinforcing bars in row one is NI. This number is 
data supplied by the user. 
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The minimum number of bars in row one or the last 
row is two. The program itself sets no limitation 
upon the maximum number of bars which may be in 
the row. This limitation is a practical one set by 
feasability of good and uniform construction. There 
is no program check as to whether the bar spacing 
supplied by the user is below the minimum allowed. 
There is also no program check as to whether the 
number of bars specified by the user may be included 
within the dimension DX with a cover distance of CV 
and bar spacing of CBX. The user must make this 
check himself. 

For reinforcing bars in row one RY is constant and 
equal to CV. RY is the distance to the barri question 
from DC . The key as to a bar being in tension or com- 
pression lies in its relationship to RPX1 . If RX (the 
distance from AD to the bar in question - see fig. 1 & 2) 
is less than RPX1 the bar is in tension. If RX is greater 
than RPX1 the bar is in compression. 

Reinforcing for rows 2- (NJ-1) are considered next. NJ 
is the number of rows of reinforcing steel and is speci- 
fied by the user. For rows 2-(NJ-l) only two bars of 
steel exist, one in column one and one in column NI. 

If the bar under consideration is in column one, the 
value of RX is CV. The critical dimension is RPY1 . 
If RY is greater than RPY1 the bar is in tension. If RY 
is less than RPY1 , the bar is in compression. 

If the bar under consideration is in column NI, the 
value of RX is constant for all bars in the column. The 
critical dimension is RPY2. If RY is less than RPY2, 
the bar is in compression. If RY exceeds RPY2, the bar 
is in tension. 

In row NJ the value of RY is constant for all reinforcing 
bars in the row. The key dimension is RPX2 . If RX 
is less than RPX2, the bar is in tension. If RX exceeds 
RPX2 the bar is in compression. 

If a bar is in tension it's area is multiplied by the ratio 
AN in computing its contribution to the transformed 
section. AN is the ratio n or ratio of flectural modulus 
of elasticity for reinforcing steel to flectural modulus of 
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elasticity for concrete . For concrete having a 
compressive strength of 3000psc @ 28 days n is 
10. If a bar is in compression, its* area is multi- 
plied by (n-1) to compute its contribution to the 
transformed area. For each bar the following com- 
putations are done in phase one. First, its trans- 
formed area is added to the sum of the transformed 
areas (SUMA). Next, the moment of its transformed 
area about AD is added to the sum of transformed area 
moments about AD (SUMAX). Finally, the moment of 
its transformed area about DC is added to the sum of 
transformed area moments about DC (SUMAY) . 

After all steel reinforcing has been processed, the 
concrete area in compression is considered. In 
addition to its affect upon SUMA, SUMAX, and SUMAY 
the moment of inertia about X and Y axis through the 
centroid of concrete areas in compression is computed 
at this time . A glance at figure 3 of appendix shows 
three possible concrete areas. The program tests for 
the presence of each type and calculates its contri- 
bution . Each of the areas present has its own centroid . 
It is about the X and Y axis of these centroids that we 
compute Moments of Inertia at this time. 

XBAR (the location of YY axis referred to AD) is com- 
puted by dividing SUMAX by SUMA. YBAR (the location 
of XX axis referred to DC) is computed by dividing 
SUMAY by SUMA. The revised moments about the XX 
and YY axis are now computed. These will differ from 
the original moments if XBAR and YBAR differ from DXE 
and DYE (assumed dimensions from AD and DC to point 
of action of axial load) . 

After locating the XX and YY axis for the column, the 
program is now able to compute the contribution of the 
steel reinforcing to the moments of inertia about the 
XX and YY axis . Once again a determination is made 
as to whether a bar is in tension or compression. 

The concrete areas in compression contribute further 
to the moments of inertia because of the eccentricity 
or distance of their centroids from the centroid of the 
transformed section. This contribution equals the area 
times the distance squared . 
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Phase three of the program covers the actual 
computation of stresses, now that the character- 
istics of the transformed section have been cal- 
culated. The stress due to axial load (P/A) and 
the stresses due to moments (MC/l) are now com- 
puted . MC/l is determined for all four column 
edges. The stress at any corner is the sum of 
P/A, MyCy/ly, and MxCx/lx. The stress for the 
bar in corner A (see fig.l) is computed. This is 
the bar having maximum steel stress. Part of the 
input supplied to the program by the user were four 
dimensions describing the location of the neutral 
axis assumed for purposes of analysis. The location 
of the neutral axis is computed in phase four based 
upon the stresses computed in phase three. If the 
computed location agrees closely with the assumed 
location, the section characteristics and stresses 
computed are close approximations of the actual 
values. If there is a large discrepancy, the section 
characteristics and stresses differ widely from the 
actual. The two examples shown in appendix 
bear this out . 

The computer will stop calculations after phase four. 
The printed output obtained at all phases is shown 
in the description of typed output records on pages 8 
to 9 . The last line given shows the computed neutral 
axis location. Should this differ widely from the 
assumed and should the user wish to perform another 
analysis, he need only push start. The computed 
neutral axis location will be used as the assumed 
neutral axis location in a new analysis. This pro- 
cedure may be continued as often as required. 
Example 2 shows that the number of iterations re- 
quired is not large even when a bad guess is made 
as to neutral axis location for our initial assumption. 



Typed Output Records 

(see appendix for symbol explanations) 



Record 1 SUMA 
SUMAX 
SUMAY 

Record 2 AC (3) 

SUMAY 
SUMAX 
SMRIY 
SMRDC 

Record 3 AC (2) 

SUMAY 
SUMAX 
SMRIY 
SMRDC 

Record 4 AC (1) 

SUMAY 
SUMAX 
SMRIY 
SMRDC 

Record 5 XBAR 
YBAR 

Record 6 SUMA 
0MYR 
0MXR 

Record 7 SMRIY 
SMRDC 

Record 8 ARKY 
ARKX 
SMRIY 
SMRDC 

Record 9 ARKY 
ARKX 
SMRIY 
SMRDC 

Record 10 ARKY 
ARKX 
SMRIY 
SMRDC 



this record typed out only if concrete area 
3 present ( see fig. 3) 



this record typed out only if concrete area 
2 present (see fig. 3) 



for concrete area L 



for concrete area 3 if present 



for concrete area 2 if present 
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Record 11 



0MIY 

0MDC 

TSADY 

TSBCY 

TSDCX 

TSABX 



Record 12 



SA 

SB 
SC 
SD 
SST 



Record 13 



DDYI 
DDXI 
DBYI 
DBXI 



typed at end of each iteration 
computed neutral axis location 
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PROGRAM DESCRIPTION -- LISTING PROGRAM 



PURPOSE: 



INPUT: 



To save time the detail Information concerning 
steel reinforcing, its effect upon the location 
of XBAR and YBAR, and its contribution to moment 
of inertia about X & Y axis is punched into cards . 
This program is to list this detail data if desired. 

Output card types 1 , 2 and 3 from the column 
analysis program. 



OUTPUT: 



RPX1, RPY1, RPX2, RPY2 on line 1. J, I, M, 
Factor, X & Y location of each bar. This is 
output from the computation of XBAR and YBAR. 
J, I, M, Factor, X & Y distances of each bar 
from the X and Y neutral axis . 
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Programming Notes: 



1. Origin of the object program is 8000. 

Address at end statement is 17442 (see Appendix C). 
Address at symbol table end is 18319 (see Appendix C). 
There is therefore little room for additions to the program. If 
changes are made, it is imperative that the compiler be altered 
to begin origin at 8000. Technical Bulletin No. 84 for the 1620 
covers the required changes . 

2 . Compilation and running assumes the availability of the automatic 
divide feature. Compilation for use on a machine without auto- 
matic divide requires origin at 8300 and will result in overlap con- 
dition in the present form of the source program. 

3. The object program does not have the subroutines compiled. These 
must be loaded at running time after the object program is loaded. 

4 . Corner A must be in tension. 
Comer C must be in compression. 

Running the Program ( Main Program) 



1 . Clear Memory 

Reset 
Insert 

26 00008 00009 

Release 

Start 

After 2 seconds , instant stop 
Reset 

2 . Place obj ect program in reader hopper . 
Press Load Key of card reader. 

3. The reader will stop with the last two cards unread. 
Hit reader start to permit loading of the last two cards . 

4. After the object deck has been loaded, the message "Enter 
Subroutines, Push Start" will be typed. 

5. Place Subroutines Deck with Automatic Divide in reader hopper. 
Press reader start and console start. See (3). This also applies 
here. The message "1620 Fortran Sub. Auto. Div. 9/30/61" 
will be typed. 
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6. After the subroutine deck has been processed and the required 
subroutines have been loaded, the message "Load Data" will be 
typed . 

7 . Place data cards in the reader hopper . 
Press Punch Start. 

Press reader start and console start. 

8 . See (3) ... This also applies here . To avoid this condition two 
blank cards may be inserted after the three data cards . When the 
reader stops all data will have been read and processed. Then 
hit non-process runout to clear read feed. 

9. The computer will perform the analysis of the column using the 
assumed neutral axis location read from the input data cards. A 
neutral axis location will be computed based upon the concrete 
stresses computed for the four corners A,B,C, and D. The com- 
puter will halt at the end of the first iteration . 

10. Should another iteration be desirable, press console start. Another 
analysis of the column will be performed using the neutral axis com- 
puted on the previous iteration as the assumed neutral axis. 

11. Step 10 may be repeated as often as desired. 

12. After the last iteration has been completed, lift remaining blank 
cards from punch feed and hit non-process runout. Remove output 
cards from stack and hold. The last two cards are blanks. 

13. To process another column,, two options are available. 

14. If data cards are placed in reader feed for one column at a time, 
place data cards for new column in reader hopper. Press punch 
start. Press reader start and console start. Steps 8-12 are 
applicable here. 

15. If data cards are placed in reader feed for more than one column at 
a time, step 8 applies only to the last column to be analyzed. To 
start analysis of next column press punch start and console start. 
Steps 9-12 are applicable here. 

Running List Program 
* 

1 . The output cards removed in step 12 of operation of the main program 
contain detail data on the contribution of the steel reinforcing during 
phases I and III. This program will list these cards if it is so desired. 

♦Control card must be inserted before data cards (see pg. 12) 



Steps 1 - 8 of the operation procedures for the main program 
apply here also. 

The computer will halt after listing the output cards from the 
first iteration. To list the cards from the succeeding iteration, 
press start. 

Step 3 should be repeated until output cards from last iteration 
have been listed. 

To list output cards from the analysis of another column, place 
cards in reader hopper, hit reader start and console start. 



- 20 - 





ft*** 





© 




- 23 - 




. i . 

c*/.c 

Plfi-C-fii 
S urn n 

So AM A X 

So *V1 A V ! 




T 




/ 

3 



- 24 - 



LO Z 

c5Ze~ ^~ 




Jo** 




/OS 


C^LC 






Sum/) 











XRs 1 



/IO 











S*Mt*AY> 




SMAtV, 





I 



Tinc- 



MO) r 






-cv 



JtuifP- 2. 





* — tv 





t>MAV = 

«Y-7 





3ZO 



SCCTIOnl 3 









*&CT 1 o rj | 














Ty#>e 












AX>H 


X 






smn IV 














swrfir- — 1 




o 

ll 

* 
o 
a 


► 


II o 
A £ 




n 
j 






4 


i 











II 
> 


-* 




— ♦ 


r 

it 


a 
a 











"5 

o 





1! 




ift 






COMPUTER y 
TECHNOLOGY 



- 32 - 



260000800009 



DIMENSION AC(3) 
3300 FORMAT (2H ) 

1001 FORMAT (F6.2,F6.2,F6.2,F6.2,F5.2,F12.2,F12.2,F12.2,I3,I3) 

K»5 FORMAT (F9.3.2H ,F9.3,2H ,F9.3.2H ,F9.3,2H ,F10.3,2H ,F10.3) 

1007F0RMAT (13, 2H ,13, 2H ,13, 2H ,F6.2,2H ,F7.2,2H ,F7.2) 

%mk FORMAT (F10»0^2M .F10.0.2H .F10.0.2H ,F10.f),2H ,F10.0) 

C READ IN COLUMN DATA 

REAP 1001 , 0*y~QX, CV, CBX, CBY, DXE, DYE 

READ l65l , DDY1, DDX1 , DBY1, DBX1, AS, P, OMX, OMY, Nl, NJ 
-READ tO0l T AH, ANMt 
3KK - 1 

_ . KR - 1 

SUMA - 0.0 

_ SUMAX ■ , 

SUMAY « 0.0 
C CALCULATE XlQQRDiNATES OF PI AND P2 
YDIF - DY - DBY1 - DDY1 
XPIF - DX ^-DBX4 - DDX1 
IF (DDY1) 12,14,12 

I 2 RP X 1 » CV ■-■ 

I^Y^-^DYI + CV*YDIF/XDIF 

14RPY1 »~CV 

RP X1 m DDX 1 + CV*XD IF/YD IF 
• 16IF (DBYl) 18,20,18 

RPY2 - OY - DBYl - CV*YDIF/XDIF 
GO TO 22 

20RPX2 - DX - DBX1 - CV*XD IF/YD IF 

RPY_2_-_0Y - CV 
22 PUNCH 1005, RPX1, RPY1, RPX2, RPY2 
C CALCULATE AREA OF STEEL AND MOMENT OF STEEL AREA (AY,AX) 
DO 697 J - 1,NJ 
IF U-i) 26*28^26 
26IF (J-NJ) 30,32,30 
JQM-- 3 
1-1 



4902262 

GO TO 35 
32M * 2 

GO TO 34 
28M - 1 

34DO 612 I - 1,NI 
35 KWIP - 4 
36V - I 

RX - CV + (V-1.0)*CBX 

V - J 

RY - CV + (V-1.0)*CBY 
C TEST IF BAR IS IN COMPRESSION OR TENSION AND SET FACTOR 

GO TO (48, 68, 58,38 ),M 
38IF (RY - RPY2) 40,42,42 
40FACTR - ANM1 



e 



GO TO 80 . - 33 - 0520 

42FACTR m AM ' 0530 

GO TO 80 0540 

48IF (RX - RPX1) 50,52,52 0550 

50FACTR = AN 0560 

GO TO 80 0570 

52FACTR « ANM1 0580 

GO TO 80 0590 

58IF (RY -RPY1) 60,62,62 0600 

nnin 60FACTR «= ANM1 0610 

np,n G0 T0 80 0620 

nnan 62FACTR « AN ' 0630 

nnfn G0 T0 80 0640 

00^0 681F (RX - RPX2) 70,72,72 0650 

0060 70FACTR = AN 0660 
0070 
0080 
0090 
0100 
0110 
0120 

nian auriMA = ouhaa t- mli.~ra-- w/ju 

Ol4o SUMAY " SUMAY + ACC*RY 0740 
0150 
0160 
0170 
0180 



GO TO 80 0670 

72FACTR = ANM1 0680 

80 GO TO (30 '), 214,214, 81 )>KWI P 0690 

81 PUNCH 1007, J, I, M, FACTR, RX, RY 0700 
ACC «= AS*FACTR 0710 

SUMA = SUMA + ACC 0720 

SUMAX = SUMAX + ACC*RX 0730 



GO TO (612,612,90,697) ,M O75O 

90 M = 4 0760 

I » Nl O77O 

GO TO 35 0780 

si »« 0^ 

0200 c CALCULATE AREA AND MOMENTS OF AREA FOR CONCRETE 0810 

02 »° PRINT 1008, SUMA, SUMAX, SUMAY 0820 

0220 SMRIY = 0.0 0830 

0230 SMRIX - 0. 0840 

0240 IF (DDY1) 102,104,102 0850 

0250 102AC(3) = 0X*DDY1 0860 

0260 SUMA - SUMA + AC (3) 0870 

0270 SUMAX = SUMAX + AC(3)*DX/2.T 0880 

0280- SUMAY = SUMAY + AC(3)*DDY1/2.0 0890 

0290 SMRIY «= AC(3)*DX**2/12.0 + SMRIY 0900 

0300 SMRIX - AC(3)*DDY1**2/12.0 + SMRIX 0910 

0310 KK - 2 0920 

0320 PRINT 1008, AC(3), SUMAY, SUMAX, SMRIY, SMRIX O93O 

0330 104IF (DBX1) 108,110,108 0940 

0340 108AC(2) ■= DBXlHOY - DDY1) 0950 

0350 SUMA = SUMA + AC(2) O96O 

0360 SUMAY = SUMAY + AC(2)*(DY + DDY1)/2.0 0970 

0370 SUMAX = SUMAX + AC(2)*(DX - DBX1/2.0) O98O 

SMRIY - SMRIY + AC (2)*(DBX1**2>/12.0 O99O 

SMRIX m SMRIX + AC(2)*(DY-0DY1)**2/12.' 1000 

KR - 2 1010 

PRINT 1008, AC(2), SUMAY, SUMAX, SMRIY, SMRIX 1020 

110SX - DX - DDX1 - DBX1 1030 

0380 SY - DY - 0DY1 - DBYl 1040 

0390 ac(1) = sx*sy/2.h 1050 

0400 SUMA - SUMA + AC(1) 1060 

0410 SUMAX « SUMAX + AC(1)*(DDX1 + 2.0*SX/3.0) 1070 

0420 SUMAY «= SUMAY + AC(1)*(DDY1 + SY/3.^>) 1080 

0430 SMRIY «= SMRIY + AC (1)*SX**2/18.n 1090 

0440 SMRIX m SMRIX + AC ( 1 )*SY**2/18.0 HOO 

0450 PRINT 1008, AC(1), SUMAY, SUMAX, SMRIY, SMRIX 1110 

0460 XBAR = SUMAX/SUMA 1120 

0470 YBAR = SUMAY/SUMA 1130 

0480 OMYR = OMY - P*(XBAR - DXE) 1140 

0490 OMXR - OMX - P*(DYE - YBAR) 1150 

0500 PRINT 1005, XBAR, YBAR 11 60 

0510 PRINT 1008, SUMA, OMYR, OMXR 1170 



C CALCULATE DISTANCE TO THE NEUTRAL AXIS OF STEEL - 34 - 
C AND OBTAIN THE FACTOR 
DO 967 J m 1,NJ 
IF (J-l) 202.206,202 
202 IF (J-NJ) 208,204,208 
20ADNAY - DY - YBAR - CV 
H - 2 
GO TO 210 
i06DNAY - YBAR -CV 

M « 1 
210D0 912 I ■ 1 ,N I 
. V . 1 
RX - CV + (V-1.0)*CBX 
DNAX - RX - XBAR 
KWIP - 1 
GO TO 36 
212DNAX - DX - XBAR - CV 
KWIP ml 
GO TO 36 
20 8M - 3 

DNAX - XBAR - CV 
KW1P m 3 
GO TO 36 
_2illV.» J__ 

RY - CV + (V-1.0)*CBY 
DRAY m RY - YBAR 
300 PUNCH 1007, J, I, M, FACTR, DNAX, DNAY 
£_ ^CALCULATE MOMENT OF INERTIA FOR STEEL ABOUT X AND Y AXIS 

ARKY « AS*(FACTR*(DNAX*DNAX) ) 
. ARKY *_AS*tEAC TR* ( DNAY*D NAY ) ) 
SMRIY - SMRIY + ARKY 
SMRIX - SMRIX + ARKX 
GO TO (912,912. 302, 967), M 

jQ2M m 4 

GO TO 212 
912rf>NTINIIF 
967C ONT I NUE 

-PRINT 1008, SMRIY, SMRIX 
C CALCULATE AD SOD FOR CONCRETE 

_^RKY - AC(1)*(DDX1 + 2.0*SX/3.0 ~ XBAR)**2 
SMRIY - SMRIY + ARKY 
_ .ARKX- ACXU*CD0Y1 + SY/3.0 - YBAR)**2 
SMRIX - SMRIX + ARKX 
PRINT 1008*. ARKY, ARKX, SMRIY, SMRIX 
GO TO (322,320) ,KK 
320ARKY - AC(3)*(DX/2.0 - XBAR)** 2 
SMRIY » SMRIY + ARKY 
ARKX - ACX3)*(YBAR-DDYl/2.0)**2 

SMRIX - SMRIX + ARKX 
PRINT 1008, ARKY, ARKX, SMRIY, SMRIX 
322G0 TO (332,330) ,KR 

3 |0 ARKY - AC(2>*(DX - XBAR - DBX1/2.0)**2 
SMRIY - SMRIY + ARKY 

ARKX - AC(2)*(DY - YBAR - (DY-DDY1 )/2.0)**2 
SMRIX - SMRIX + ARKX 
PRINT 1008, ARKY, ARKX, SMRIY, SMRIX 
C CALCULATE CONCRETE AND STEEL STRESS 
132POA - (-P/SUMA) 
OMIY - OMYR/SMRIY 
0M1X - OMXR/SMRIX 
TSADY - OM I Y*XBAR 
TSBCY - (-OMIY*(DX - XBAR)) 
TSDCX - ( -OM I X*YBAR ) 
TSABX - OMIX*(DY - YBAR) 

PRINT 1005, OMIY, OMIX, TSADY, TSBCY, TSDCX, TSABX 



1180 = n)h + t«apy + ii>At.x - 35 - ' a 3° 

119C SB = POA + TSBCY + TSABX ' 1840 

1200 SC = POA + TSBCY + TSDCX 1850 

121C SD = POA + TSADY + TSDCX i860 

1220 SST = ((0MIY*(XBAP.-CV))*AS + POA + AS*(OMIX*(DY-YRAR-CV) ) )*10 .0 1870 

1230 PRINT 1005, SA, SB, SC, SD, SST 188C 

1240 IF (SB) 501,503,503 1890 

1250 501DBY1 = 0.0 1900 

1260 DOUBT = DX*SB/(f>A-SB) 1910 

1270 L = 1 1920 

1280 GO TO 2000 193C 

1290 3001 DBX1 = DOUBT 1940 

1300 GO TO 505 1950 

1310 503DBXI =0." I960 

1320 DOUBT = DY*Sf'./(f»B-SC) 1970 

1330 L = 2 1980 

1340 GO TO 2000 1990 

1350 3002 DBY1 = DOUBT 2000 

1360 505IF (SD) 507, 504, 50 q 2010 

1370 507DDX1 = 0.0 2020 

1380 DOUBT = DY*SD/(SA-SD) 2030 

1390 L = 3 2040 

1400 GO TO 2000 2050 

1410 3003 DDY1 = DOUBT 2060 

1420 GO TO 511 2070 

1430 509DDY1 =0. 2080 

1440 DOUBT = DY*S[)/(SD-SC) 2090 

1450 L = 4 2100 

1460 GO TO 2000 2110 

1470 3004 DDX1 = DOUBT 2120 

1480 511 PRINT 1005, DDY1, DDX1, DBY1, DBX1 2130 

1490 PAUSE 2140 

1500 PRINT 3300 2150 

1510 GO TO 3 2160 

1520 2000 IF (DOUBT) 2001,2002,2002 2170 

1530 2001 DOUBT = D0UBT*(-1. >) 2180 

1540 2002 GO TO (3001 ,3002,3003,3004) ,L 2190 

1550 END 2200 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 
1670 
1 680 
1 690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 



- 36 



260000800009 

ENTER SOURCE PROGRAM, PUSH START 

T7358 2002 GO TO (3001, 3002, 3003, 3004), L 

T7kk2 END 

PRO G SW 1 OttfOR SYMBOL TABU, PUSH START 

T9999 SIN 

T-99B9 SJNF 

T9979 COS 

T99&9 COSF 

T9959 ATAN 

T99*»9 AT'ANF 

T9939 EXP 

T9929 EXPF 

19919 LOG 

19909 LQGF 

T9899 SORT 

19889 -SQRTF 

T9879 AC T9859 

I98&9 1300 

T°°39 3300 

J. 29 Toot _ 
T98t9 Tool 
T980 9 T005 
T9799 t005 
T97fl9 1007 
T9779 T007 
T9769 T008 
T9759 T008 
T97 fc 9 OY - 
T9739 DX 
J972gxv 
T9719 CBX 
T9709 CBY 
T9699 DXE 
- T9689 -DYE 
T9679 DDYl 
49669 WW J 
T9659 0BY1 
T9649 DBX1 
T9639 AS 
T9629 * 
T9619 OMX 
T9609 OMY 
T9599 «» 
T9589 WJ 
T' 7 9 AN 
f9*(>9 AI4M1 
T9559 5003 
T95^9 KK 
T9539 50 OT 
19529 KR 
T9519 SUMA 
Tqqoq 0ononnnn33 



19^79 
T94fa c ; 

T9459 

T9449 

Jo' oq 

Ti, -9 
T9409 
T9399 
T9389 
T9379 
T9369 
T9359 
T9349 
T9339 
T9329 
T9319 
T9309 
T9299 
T9289 
T9279 
T9269 
T9259 
T9249 
T923S 
T9229 
T9219 
T9209 
T9199 
T9189 
jo-t 9 

TS >9 
T9159 
T9149 
T9139 
T9129 
T9119 
T9109 
T9099 
T9089 
T9079 
T9069 
T9C59 
T9049 
T9039 
T9029 
T9019 
T9009 
T8999 
T8989 
T8979 
T8969 
T8959 
T8949 
T8939 
T8929 

T£ 9 
T8909 

T8899 



SUMAX 
SUMAY 
YD1F 

000 

XDIF 
CO 12 

no 1^ 

RPX1 

RPY1 

0016 

0018 

CO 20 

RPX2 

RPY2 

001 

CO 22 

0697 

J 

0026 
0028 
0030 
0032 

0003 
I 

0035 
0002 
0034 
0612 
KW1P 
0005 
5036 
V 

RX 

Toooo 00001 

RY 

0048 

0068 

0058 

0038 

0040 

0042 

FACTR 

0080 

0050 

0052 

0060 

0062 

0070 

0072 

0300 

0214 

0081 

ACC 

0090 

SMRIY 

SMR I X 

0102 

0104 

200000000 1 
T 200 000 CO 2 
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T8889 0108 



I vsu/ .< v. ■ I V 

T8869 SX " 38 



T6859 SY 
T8849 lOCOOOOOCl 
18839 18OOOO0002 



T ^829 XBAR 
.o8l9 YBAR 
T8809 OMYR 
T8799 OHXR 
T8785 S967 
T8779 0202 
18769 5206 
T8759 0208 
187^9 5204 
18739 DNAY 
T8729 6210 
T8719 0912 
I8709 DNAX 
T8699 0212 
18689 ARKY 
186 7 9 502 
T0669 "ARKX ■ 
T8659 0302 
18649 0322 
T8639 O320 
T8629 §332 
18619 0330 
T8605 0©3 
T8599 POA 
18589 OWY 
~9579 OMIX 
18569 TSAOY 
T8559 TSBCY 
185*9 TSOCX 
T8539 TSABX 
T8529 SA 
T8519 SB 

T8509 sc 

T8499 SD 

imts SST 
T8479 006 
T8469 T000000052 
T8459 05O1 
T8W 5503 
18439 DOUBT 
f 8429 1 
184 19 2000 
18409 5001 
T8399 0505 
T8389 5002 
T8379 0507 



18369 C509 
18359 3003 
T8349 0511 

T8339 3004 
T8329 Z001 
18319 2002 

SW 1 OFF TO IGNORE SUBROUTINES, PUSH START 
PROCESSING COMPLETE 



© 
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ENTER 

08000 

08036 

08096 

08156 

08180 

082C4 

08228 

08324 

08336 

08420 

08540 

08564 

08588 

08612 

08624 

08708 

08828 

08840 

08848 

08876 

Y ) 

09014 

09116 

NAY ) 

09254 

09340 



1 



SOURCE PROGRAM, PUSH START 
READ 1003, Nl , NJ 
READ 1005, RPX1, RPY1, RPX2, RPY2 
PRINT 1005, RPX1, RPY1, RPX2, RPY2 
PRINT 1001 
PRINT 1002 
PRINT 1001 

KUK = 2*NI + 2*(NJ-2) 
DO 77 N = 1.KUK 

READ 1003, J, I, H, FACTR, RX, RY 
PRINT 1003, J, I, M, FACTR, RX, RY 
PRINT 1001 
PRINT 1004 
PRINT 1001 
DO 88 N = 1,KUK 

READ 1003, J, I, M, FACTR, DNAX, DNAY 
PRINT 1003, J, I , M, FACTR, DNAX, DNAY 
PAUSE 
GO TO 1 

FORMAT (2H ) 

FORMAT (5H J ,5H I ,5H M .8HFACT0R 



77 



88 



1001 
1002 



1003 
1004 



,9H RX 



,7H R 



FORMAT (I3,2H 
FORMAT (5H J 



.I3.2H ,I3,2H ,F6,2,2H .F7.2.2H t F7.2) 

,5H I ,5H M ,8HFACT0R ,9H DNAX ,7H D 



1005 FORMAT (F9.3.2H ,F9.3,2H ,F9.3,2H ,F9.3,2H ,F9.3) 

END 



PROG SV.' 1 ONFOR SYMBOL TABLE, PUSH START 

T9999 SIM 

1998^ SIMF 

T9973 COS 

T9969 COSF 

T9959 ATAN 

T9949 ATANF 

T9939 EXP 

T9929 EXPF 

T9919 LOG 

T9909 LOGF 

T9899 SORT 

T9889 SORTF 

T9879 TC03 

TS869 T0G3 

T9859 Nl 

T9849 f!J 

T9839 Ocoi 

Tq82<i To OS 



T9819 T005 - 40 

19809 RPX1 

9799 RPY1 
RPX2 
*9779 RPY2 
TOO I 

9759 1001 
19749 1002 
T9739 TOO 2 
T9729 KUK 
T97t9 0002 
19709-1100 
T9699 002 
19689 0077 
f 9679 N 
19669 J 
T9659 I 
19649 M 
T9639 FACTR 
19629 RX 
T9619 RY 
I96Q3 iDQk 
T9599 TOO** 
T9589 Onflfl 
T9579 DNAX 
19569 DNAY 



SW I OFF TO IGNORE SUBROUTINES, PUSH START 
PROCFSSillG.-CQMP.LETE 
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Appe^ot/. fc- TEST * I 



7200 7200 350 
1020 000 1000 
1000 900 



600 600 3600 3600 

000 156 113200000 *f536000000 3070800000 12 12 
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CTTER SUBROUT \W*> , PUSH START 

1620 FORTRAN SUBR. AUTO DIV 9/30/61 

LOAD DATA 



EKflJH puts- I 



1ST 



+550. 
+734. 
+1864. 
+42.843 
+32*9. 
+t3g6l7. 

♦^39*1 
+IE.5B3 

+12.^91 



2MD 



^55TT 
+899V 
+I45BT 
+41.821 
+3021. 
+1303112. 
♦Mot?; 

♦So484. 
+17.355 
^1 830. 501 

+ir.r 



+24394. 
+49719. 
+139229. 

+32252003. 



+23280. 
+28139. 
+79359. 
+24.420 
+22961136. 
+852 102. 

+274138. 
+28.203 
+527.564 -1503.090 
. +.000 +18.705 



+317260. 
+854323. 



+6367. 
+284350. 



+1402201. +869406. 
+1436596. +1143544. 
+684.768 -466.008 -688.744 
-352.313 +25720.226 
+.000 



+1341.910 



+53457^ 
+39 182. 
+68509T 
+2 2.670 
+24T17723. 



+74961 
"~ +1720-. 
+242636. 
+29i587 
+580.936 
+;oo<r 



+24564. 
+55876. 
+126383. 

+30271212. 



+388545. 
+811586. 



+11695. 
+147556. 



+1359T80. +766816. 
+1389664. +1009453. 
+725*820 -523.744 -679.842 
-1578.180 -328.615 +28068.623 
+19.372 +.000 



H1479.274 



3RD 



"+553T +21497. +24564 . 
.+788. +28887. +51901. +340847. 
+1560. +66161. +123905. +772874. 
+42.105 +22.483 
+2942. +23796303. +30058965. 
+758265. 

+8389; +1318831. +766655. 
+228128. +1348245. +994783. 
+30.216 +743.159 -527.628 -679.364 
+583,923 -1591.670 -320.882 +28474.175 
+.000 +19.324 +.000 



+7895. 
+152598. 



♦29413. 
+17.649 
+1854.711 
+10,619 



+1496.228 



o 

- 43 - 

7200 72CO 350 600 600 3600 36OO 

00) 1000 000 1000 156 113200000 4536000000 307O8OOOOO 12 12 
1000 900 
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CJZ. Con- 

TGS T NO x 



- 45 - 



ENTER SUBROUTINES, PUSH START 

1620 FORTRAN SUBR. AUTO DIV 9/30/bt 

LOAD RATA 



4908000 



TFT 



♦652. 
♦720. 

♦47.765 

+•$5191. 

+J757JV 
+266377. 
+15v26fl 
♦1503. 



♦23098. 
+50132. 
+95Q60. 
♦29.302 
+173JJ9533. 

+32293. 

+26.316 



+24212. 
+71338. 
+154954. 

♦37778815. 

+873167. 
+113954%. 
•728.903 



+6000. 
+287216. 



(T y A *» J» i- C- 



+311040. 
+850176. 



+12.527 



+3.500 



503.592 +404.867 -1489.907 
♦14.164 +..QPQ +15384 




+1403273. 
+1435566. 
-369.821 
-391.182 
♦.000 



-771.142 
+23140.053 



+1123.631 



♦23383V +24460. 

+32414-. +61912. +462349. +19706. 
+75673. +134057. +895220. +165259. 
♦23T473 
+24*8708 3 . +31 180010. 
♦751065. 

♦42342V +1443369. +793413. 
+275397. +1*76739. +106881© . 
+22kl72 +686.721 -502.296 -684.783 +1415.644 
+562.209 -1538.218 -349.200 +26791.033 
+^Q00__ .+13.271 +.000 



♦121 
♦54744. 
+29917. 
+17.525 
+1844.715 
. +10.763 



+13437i 
+29722. 
+62221. 
+22.5; 
+21 

+13380, 
♦236100. 



+24564. 
♦54524. 
♦124953. 

v+JOl 19178. 



+372319. 
+794898. 



+10290. 
+145706. 



+1341729. 
+1371647. 



4 +734.165 
-1586.061 
+J9.349 



+763726. 
+999826. 



-527.654 
-324.241 
+.000 



-678.894 

+28301.151 



+1490.062 



+653;. 
+774. 
+1507. 
+42.150 
+2936. 
+1263286. 
+51593. 
+29317. 
+17.665 
* 855. 947 
+10.602 



+23497. +24564. 

+28735. +51396. +334785. +7481. 
+66012* +123776. +769069. +154464. 
+22.482 
+237i5M5- +30057879. 
+76OI50. 

+7589. +1314879. +767740. 
+226621. +1344197. +994362. 
+30.228 +744.598 -527.293 -679.599 +1496.837 
+584.056 -1592.380 -320.489 +28496.534 
+.000 +19.321 +.000 



J 


I 




FAC TOR 


4-J 


+1 


4-| 


+tc .00 


+ 1 


+7 
+* 


^ 1 


+ 10 00 


4-1 


+3 


+f 


x() 00 




+Ji 


4.* 


+9 00 


4.* 


4.C 


4.| 


+9.00 


^ f 


+6 


j_l 
' V 


+<J )0 


■ | 


4.7 


4-| 


+9 . 00 


+ t 


+fi 


4.] 


+9 00 




4-4 


4-1 


+9 00 


^ 1 


+lfi 

T|U 


4.1 


+9.00 


' 1 


'II 


4>| 


+9.00 


' I 


4-17 

^ W ML 


4.I 


+9 00 


+2 


4.| 


4.2 


+10.00 




+1? 


+4 


4.0 no 


+a 


+| 


+* 


+10.00 


+a 


+12 


+4 


+9.0 1 






+3 


+10.00 


+4 


+12 


+4 


+9,00 




4vf 


+3 


+10.00 




+12 


+4 


+9 00 


+6 


+| 


+3 


+10.00 


+6 


+12 


+4 


+9 00 


■ / 


+f 


+a 


+10 . >o 


+7 


+12 


+4 


+5.00 


+8 


+1 


+3 


+10.00 


+8 


+12 


+4 


+9.00 


+9 


+1 


+3 


+10.00 


+9 


+12 


+4 


+S.00 


+10 


4.) 


+a 


+10.00 




+17 


+4 


+9 . 00 




4.1 


4.3 


+10 .00 


• it 


+17 


+4 


+9 00 


+12 


+ 1 


+2 


+10.00 


+12 


+2 


+2 


+10.00 


+12 


+3 


+2 


+10.00 


+12 


+4 


+2 


+10.0 5 


+12 


+5 


+2 


+1C.00 


+12 


+6 


+2 


+10.00 


+12 


+7 


+2 


+10.00 


+12 


+8 


+2 


+10.00 


+12 


+9 


+2 


+10.00 


+12 


+10 


+2 


+10.00 


+12 


+11 


+2 


+9. ;1 


+12 


+12 


+2 


+9.00 


J 


1 


1-1 


FACTOR 


+1 


+1 


+1 


+10.00 


+1 


+2 


+1 


+10.00 


+1 


+3 


+1 


+9.00 


+1 


+4 


+1 


+9.00 


+1 


+q 


+1 


+o„ n o 



+59.472 

RX 

+3.50 

+9.50 
+15.50 
+21.50 
+27.50 
+33.50 
+39.50 
+45.50 
+51.50 
+57.50 
+63.50 
+69.50 

+3.50 
+69.50 

+3.50 
+69.50 

+3.50 
+69.50 

+3.50 
+69.50 

+3.50 
+62.50 

+3.5© 
+69.50 

+3.50 
+69.50 

+3.50 
+69.50 

+3.50 
+69.50 

+3.50 
+69.50 

+3.50 

+9.50 
+15.50 
+21.50 
+27.50 
+33.50 
+39.50 
+45.50 
+51.50 
+57.50 
+63.50 
+69.50 

DNAX 

-44.26 
-38.26 
-32.26 
-26.26 
-20.26 



+68.500 

RY 

+3.50 

+3.50 

+3.50 

+3.50 

♦3.50 

+3.50 

+3.50 

+3.50 

+3.50 

+3.50 

+3.50 

+3.50 

+9.50 

+9.50 
+15.50 
+15.50 
+21.50 
+21.50 
+27.50 
+27.50 
+33.50 
+33.50 
+39.50 
+39.50 
+45.50 
+45.50 
+51.50 
+51.50 
+57.50 
+57.50 
+63.50 
+63.50 
+69.50 
+59.50 
+69.50 
+69.50 
+69.50 
+69.50 
+69.50 
+69.50 
+69.50 
+69.50 
+69.50 
+69.50 

DNAY 

+25.80 
+25.80 
+25.80 
+25.80 
+25.. 80 



1ST ITeR*7l»»J 



+1 


+6 


+1 


+1 


+7 


+1 


+! 


+8 


+1 


+1 


+9 


+1 


+1 


+10 


+1 


+1 


+11 


+1 


+1 


+12 


+1 


+2 


+13 


+3 


+2 


+13 


+4 


+3 


+13 


+3 


+3 


+13 


+4 


+4 


+13 


+3 


+4 


+13 


+4 


+5 


+13 


+3 


+5 


+13 


+4 


+6 


+13 


+3 


+6 


+13 


+4 


+7 


+13 


+3 


+7 


+13 


+4 


+8 


+13 


+3 


+8 


+13 


+4 


+9 


+13 


+3 


+9 


+13 


+4 


+10 


+13 


+3 


+10 


+13 




+11 


+13 


+3 


+11 


+13 


+4 


+12 


+1 


+2 


+12 


+2 


+2 


+12 


+3 


+2 


+12 


+4 


+2 


+12 


+5 


+2 


+12 


+6 


+2 


+12 


3 


+2 


+12 




+2 


+12 


+9 


+2 


+12 


+10 


+2 


+12 


+11 


+2 


+12 


+12 


+2 



+9.')0 -14.26 

+9.00 -8.26 

+9.0'1 -2.26 

+9.00 +3.73 

+9.™ +9.73 

+9.00 +15.73 

+9.00 +21.73 

+10.00 +44.26 

+9.00 +20.73 

+10.00 +44.26 

+9.00 +20.73 

+10.00 +44.26 

+9.00 +20.73 

+10.00 +44.26 

+9.00 +20.73 

+10.00 +44.26 

+9.00 +20.73 

+10.00 +44.26 

+9.00 +20.73 

+10.00 +44.26 

+9.00 +20.73 

+10.00 +44.26 

+9.00 +20.73 

+10.0 > +44.26 

+9.00 +20.73 

+10.00 +44.26 

+9.00 +20.73 

+10.00 -44.26 

+10.00 -38.26 

+10.00 -32.26 

+10.00 -26.26 

+10.00 -20.26 

+10.00 -14.26 

+10.00 -8.26 

+10.00 -2.26 

+10.00 +3.73 

+10.00 +9.73 

+9.00 +15.73 

+9.00 +21.73 



+25.30 
+25.80 
+25.80 
+25.80 
+25.80 
+25.80 
+25.80 
-19.80 
-1S.80 
-13.80 
-13.80 
-7.80 
-7.80 
-1.80 
-1.80 
+4.19 
+4.19 
+10.19 
+10.19 
+16.19 
+16.19 
+22.19 
+22.19 
+28.19 
+28.19 
+34.19 
+34.19 
+39.19 
+39.19 
+39.19 
+39.19 
+39.19 
+39.19 
+39.19 
+39.19 
+39.19 
+39.19 
+39.19 
+39.19 
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Symbol Table List 



AC Concrete areas (see fig. 3) 

DY Column dimension in Y Direction (see fig.l) 

DX Column dimension in X Direction (see fig.l) 

CV Distance to <£ of reinforcing steel from the column edge . 

CBX Bar spacing in X direction 

CBY Bar spacing in Y direction 

DXE Assumed distance to YY axis from column edge AD 

DYE Assumed distance to XX axis from column edge DC 

DDY1 Assumed distance from D along AD to neutral axis (see fig.l) 

DDX1 Assumed distance from D along DC to neutral axis 

DBY1 Assumed distance from B along BC to neutral axis (see fig.l) 

DBX1 Assumed distance from B along AB to neutral axis 

AS Area of steel reinforcing bar 

NI Number of reinforcing bars in X direction 

NJ Number of reinforcing bars in Y direction 

AN Ratio of flexural modulus of elasticity for reinforcing steel to 

flexural modulus of elasticity for concrete 
ANM1 AN - 1.0 

P Axial Load 

0MX Moment about axis XX 

0MY Moment about axis YY 

KK Switch (2 if area 3 is present - see fig. 3) 

KR Switch (2 if area 2 is present - see fig. 3) 

SUMA Sum of transformed areas 
SUMAX Sum of transformed area moments about AD 
SUMAY Sum of transformed area moments about DC 
YDIF Distance in Y direction between points where neutral axis is 

assumed to cross the column edges 
XDI F Distance in X direction between points where neutral axis is 

assumed to cross the column edges 
RPX1 Distance from AD in X direction to where neutral axis is assumed 

to cross steel row at point PI (see fig. 2) 
RPY1 Distance from DC in Y direction to where neutral axis is assumed 

to cross steel row at point PI (see fig. 2) 
RPX2 Same as RPX1 for point P2 

RPY2 Same as RPY1 for point P2 

J Counter for number of rows 

M Switch 1 for row 1 

2 for row NJ 

3 for column 1, rows 2 - (NJ-1) 

4 for column NI, rows 2- (NJ-1) 
V Temporary work area 

RX Distance in X direction from AD to bar (J,I) 
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RY 


Distance in Y direction from DC to bar (J, I) 


FACTR 


Factor relating steel area to transformed steel area 


ACC 


Transformed steel area 


I 


Counter for number of columns 


SMRIY 


Moment of inertia about XX axis for the transformed section 


SMRIX 


Moment of inertia about YY axis for the transformed section 


SX 


X dimension of triangle area 1 


SY 


Y dimension of triangle area 1 


XBAR 


Distance to axis YY from AD 


YBAR 


Distance to axis XX from DC 


0MYR 


Revised moment about YY 


0MYR 


Revised moment about XX 


DNAY 


Distance in Y direction between axis XX and centroid of steel 




reinforcing bar. 


DNAX 


Distance in X direction between axis YY and centroid of steel 




reinforcing bar. 


ARKY 


Work area for A (Dx) for concrete areas 


ARKX 


Work area for A (Dy)^ for concrete areas 


?£h 


Concrete stress due to axial load 


OMIY 


Y Moment/Y Moment of Inertia 


OMDC 


X Moment/X Moment of Inertia 


TSADY 


Stress along AD due to Y Moment 


TSBCY 


Stress along BC due to Y Moment 


TSABX 


Stress Along AB due to X Moment 


TSDCX 


Stress along DC due to X Moment 


SA 


Combined concrete stress at A (see fig.l) 


SB 


Combined concrete stress at B ( 


SG 


Combined concrete stress at C 


SD 


Combined concrete stress at D 


SST 


Steel stress for bar under maximum tensile stress ( in corner 




D - see fig.l) 


KWIP 


Switch - used to control return to phase III 


DOUBT 


Work area used in computing neutral axis location 


L 


Switch - used to control return from subroutine in phase IV 


KUK 


Counter used in listing program 



o 
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